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Abstract. - We argue that the pseudogap phase may be an attribute of the non-BCS pairing 
of quantum-critical, diffusive fermions near the antiferromagnetic quantum critical point. We 
derive and solve a set of three coupled Eliashberg-type equations for spin-mediated pairing and 
show that in some T range below the pairing instability, there is no feedback from superconduc- 
tivity on fermionic excitations, and fermions remain diffusive despite of the pairing. We conject 
that in this regime, fluctuations of the pairing gap destroy the superconducting condensate but 
preserve the leading edge gap in the fermionic spectral function. 



The pseudogap behavior in underdoped cuprates is one of the most unusual and exciting 
features of condensed matter physics. By all experimental accounts, below optimal doping 
superconducting-like behavior of fermionic observables sets in at a temperature which increases 
with underdoping. This temperature correlates with the value of the superconducting gap at 
T = 0, but does not correlate with the transition temperature itself. The latter decreases with 
underdoping and eventually vanishes jj] . 

In this paper we argue that the pseudogap behavior may be a part of new, non-Fermi- 
liquid physics associated with the pairing of fermions in the quantum-critical regime near the 
antiferromagnetic quantum critical point. The onset of the pairing in this regime has been 
recently studied by Finkel'stein and two of us (ACF) 0]. ACF demonstrated that the onset 
temperature T — Ti ns is determined by a competition between fermionic incoherence and 
the absence of a gap for spin excitations which mediate superconductivity j|. Due to this 
competition, Xi ns tends to a finite value when the spin correlation length, £, diverges. 

An obvious question posed by this result is whether the pairing instability at Tj ns leads 
to true superconductivity, or only gives rise to a formation of spin singlets which still behave 
incoherently and hence do not superconduct. If this was the case, the phase right below 
Tins would be a pseudogap phase, while the actual superconductivity would emerge only at a 
smaller temperature. 

This paper is a first step in addressing this issue. We report the results of our analysis of 
the system behavior within the Eliashberg-type approach [Q . This approach is mean-field like 
in the sense that it does not include the feedback effect from the pairing modes (e.g., phase 
fluctuations) on fermions and spin excitations. However, it does include nontrivial physics 
associated with non BCS pairing in the quantum-critical regime. We show that pairing of 
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quantum critical fermions causes a reduction of the superfluid stiffness, D s oc n s /m* . At 
low T, the reduction is predominantly due to an enhancement of the quasiparticle mass, to*. 
However, in some temperature interval between Ti ns and To < Tj„ s we found an additional 
reduction of the superfluid stiffness. We argue that the latter effect is an indication of a new 
physics associated with a quantum-critical, non-BCS pairing near the magnetic instability, 
and reflects the reduction of the superfluid density n s due to the fact that fermions remain 
diffusive at the lowest u> despite the formation of pairs. We conject that phase fluctuations, 
acting on top of the Eliashberg solution, likely destroy fermionic coherence at T ~ To, which 
then becomes a true phase transition temperature where a continuous U(l) symmetry breaks 
down. The pairing gap, defined as the scale below which the quasiparticle spectral weight is 
reduced, however, remains finite at T and disappears only at T ins . 

We now proceed with the calculations. Our theoretical treatment is based on the spin- 
fermion model which describes the interaction between low-energy fermions and their collective 
spin degrees of freedom peaked at Q = (7r,7r) The model has two energy scales: the 

effective spin-fermion interaction, g = g 2 xo, where g is the spin-fermion coupling, and xo is 
the overall factor in the static spin susceptibility \q = Xo/{£,~ 2 + (q — Q) 2 ), and a typical 
fermionic energy, vp£ , where v-p is the Fermi velocity. The dimensionless ratio of the 
two, A = 3g/(47rz)F^ _1 ), determines the relative strength of the spin-fermion coupling. The 
quantum-critical behavior obviously corresponds to the strong coupling limit, A > 1 [0]. 

For superconductors with phonon mediated pairing, a recipe to study the system behavior 
at strong coupling is the Eliashberg theory. It is justified by Migdal theorem which states 
that the corrections to the electron-phonon vertex Sg/g and vJ^dH^k, uj) / dk, although increase 
with the dimensionless coupling A, still scale as X(v s /vp) where v s and vf v s are the sound 
velocity and the Fermi velocity. As v s /vf ~ I0 -4 , Sg/g and Vp dT,(k,oS)/dk can be safely 
neglected for all reasonable A. One then has to include only E(u;) ||]. In our case, Migdal 
theorem is inapplicable as spin fluctuations are collective modes of fermions and hence the 
spin velocity (the analog of v s ) is of the same order a,s vf- This implies that the Migdal 
parameter Xv s /vp is large in the strong coupling regime. It turns out, however, that in this 
limit one again can neglect Sg/g and Vp 1 dT,{k,uj)/dk^ this time because the dynamics of 
the collective mode is completely modified by low-energy fermions and becomes diffusive at 
energies relevant to the pairing problem. [For electron-phonon interaction this happens only 
for very small frequencies ~ (v s /vf) 2 T c <C T c .] Specifically, the strong coupling solution of the 
spin-fermion problem in the normal state yields dS/dui oc A while Sg/g and Vp 1 dT,(k,uj)/dk 
scale as (l/7V)logA <C A, where N(= 8) is the is the number of hot spots in the Brillouin 
zone @j. Furthermore, below T,- ns , we found that Sg/g even becomes independent on A. 
Physically, the irrelevance of vertex corrections is due to the fact that the theory possesses 
no SDW precursors, and hence spin excitations are not near Goldstone bosons ||. Below we 
formally treat TV as a large parameter and perform computations at N — > oo. In this limit, 
vertex corrections can be totally neglected. Also, following ACF, we neglect the momentum 
variation of Ek(w) and of the anomalous vertex $k(w) along the Fermi surface, i.e., replace 
them by Ek ha (w) and $k hs (w) (respecting the d x i_ y 2 symmetry $k h3 +Q = — $k hs )- The last 
approximation is justified if g <C vpkp which we assume to hold. Finally, we verified that 
both above and below Tj ns , there is no universal thermal contribution to £ from low-energy 
fermions, i.e., £(T) ~ £(T = 0). At T > g, the physics is dominated by classical fluctuations, 
and £ -1 eventually becomes strongly T dependent [Q. 

We next discuss the structure of the Eliashberg equations for spin-mediated superconduc- 
tivity. In conventional superconductors, the Eliashberg theory involves two coupled equations 
for the fermionic self-energy E, and the pairing vertex $ jjj. Modifications of the phonon 
propagator due to fermionic pairing are small and can be neglected. For spin-mediated pair- 
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ing, the situation is different as the spin dynamics is made out of low-energy fermions and 
thus is sensitive to the opening of the fermionic gap ^,|). As a result, the Eliashberg theory 
has to be generalized to include the equation for the dynamical spin susceptibility. 

The mutual feedbacks between fermions and spin fluctuations have been earlier studied 
numerically within the FLEX approximation [^], [|. Our results agree with these studies, 
but go beyond them in the understanding of the new physics. The modification of the spin 
propagator due to fermionic pairing has also been considered in the context of marginal Fermi 
liquid phenomenology |l(J, but the feedback effect on fermions has not been studied. 

The set of Eliashberg equations for the spin-fermion model is obtained in a standard way by 
evaluating diagrammatically the pairing vertex, fermionic self-energy and the spin polarization 
operator, and integrating over the fermionic energy. The results have been earlier presented 
in the two limiting cases: for infinitesimally small $(w) Q, and for <f>(w) = const A 
straightforward extension to a general case yields ($(w m ) = $ m , etc.), 



1/2 



1/2 



E m = uj m + ^— V f" ( ^= ) , (2) 

n m = ttT ± \ (3) 



The dynamical spin susceptibility at antiferromagnetic momentum Q and the fermionic Green's 
function follow as XQ(w m ) _1 oc u si - U rn , and G k (ui m ) = -(i£ m + ek)/($m + + e 2 .), re- 
spectively, where ek = vj? • (k — k^?). For further convenience we included the bare uj n term 
in the fermionic propagator into E n and introduced <D = 9g/(2irN) and w s f = w/(4A 2 ) cx 
instead of A and g (the 1/N factor in Q is eliminated after appropriate rescaling g — > gJV and 

UF — > VfN). 

Analyzing S and %q in the normal state, we find that u s f separates the Fermi liquid 
behavior at u> < u> a i from the quantum-critical behavior which holds between uj s f and Q 3> t-u s f. 
Which of the two energies determines the superconducting properties of the system? The 
onset of pairing was determined by ACF by linearizing with respect to $. They found that 
$ emerges at a temperature Ti ns which depends weakly on uj s f and saturates at Tj ns w 0.17(D 
for £ — > oo. Thus, the onset of pairing is produced by quantum-critical fermions with 
w>w» w s f. We consider the solution of the full set of nonlinear equations ([l|-|3j) we analyze 
whether there are any new crossover scales below T ins . 

At T — 0, the solution of the Eliashberg set is in many respects similar to the one for 
$(w) = const H. The pairing gap A (defined in the usual way as A(w) = $(w)w/E(aj) and 
A(u> = A) = uj) scales with T lns up to log A corrections which, however, remain numerically 
small even for A = 20. The presence of the gap eliminates spin-fermion scattering at energies 
0(A) and restores fermionic and bosonic coherence. In particular, spin excitations, which 
were purely relaxational in the normal state become propagating below 2A with dispersion 
fi q = (A 2 + uj?(q - Q) 2 ) 1/2 , where A s ~ (w s/ ) 1 '' 2 oc T lns /X. This effect changes the 
dynamical exponent to z = 1 and yields the resonance in the inelastic neutron scattering at 
lo = A s . The spin resonance in turn affects fermionic excitations: both, ImE(w) and Im$(w) 
vanish at low u) due to a lack of a phase space for single particle decay and jump to finite 
values at uj = A + A s . This gives rise to a dip in the fermionic spectral function P,pT|. 
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Fig. 1 - The real part of the pairing vertex 3>(uj) and /mE(ai) for A = 4 and different temperatures 
T < A s , T ~ A aj A s < T < T and T < T < T ins . We associate T ~ A s and T ~ T with the 
onset of impurity-like behavior due to thermal fluctuations, and with the onset of the reduction of 
the superfluid density, respectively. 



These analytical results fully agree with our numerical solution (the lowest T results in 
Figs |l| and ||). We found 2A/T ms w 4 for A 3> 1 with downturn deviation at A ~ 1. To a 
reasonable accuracy, A s ~ 0.35 uj/\. The spectral function has a peak-dip-hump structure, 
and the peak-dip distance exactly equals A s . The numbers which we obtain are also consistent 
with experiment. We use data for Bi2212 for comparisons. Near optimal doping, which we 
identify with A ~ 1, we obtain A a /u> ~ 0.25 — 0.3 and A/cD w 0.2. The value of Q can 
be extracted from the photoemission data as a frequency where nonlinear corrections from 
Re£(ui) to the quasiparticle dispersion become irrelevant. This yields |L2| Q ~ 150— 160meV. 
Using this a), we obtain A s ~ 38 — 48meV and A ~ 30 — 32meV which are in reasonable 
agreement with the neutron scattering |]i~3| , photoemission [Il4 | and tunneling Jig] data. 

We now turn to T > 0. A simple analysis of Eqs. (0p) shows that classical thermal 
spin fluctuations (the ones with zero Matsubara frequency) critically depend on u> S f and are 
therefore relevant already at low T. These fluctuations account for scattering with zero energy 
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Fig. 2 - The dynamical spin susceptibility x"(Qi w ) an d the photoemission intensity I(fJ) = 
A(ui)riF{(jJ) for A = 4 at different T. The insets show the the superfluid stiffness D S (T) and the NMR 
relaxation rate 1/TiT vs T/Ti ns . The actual behavior of D S (T) and 1/TiT at low T is smoother due 
to contributions from nodal regions. The arrow indicates the position of the dip in I(u>). 
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transfer and therefore act for spin-mediated d-wave paring in the same way as nonmagnetic, 
elastic impurities in s-wave superconductors. We then use the same strategy as for the impu- 
rity problem |l6| and introduce $ m and E m via <I> m = $> m Vrm E m = S m 7y m , where 

ttTA 

T)m = l-\ j 

v /($ m ) 2 + (s m ) 2 

Substituting 4> m and E m into the Eliashbcrg set we obtain after some algebra that the 
quantum contributions to the self energy, $ TO and E m , obey the same Eqs. (Q), (^) but 
without zero frequency (m — n) contributions. This implies that the effects of thermal spin 
fluctuations are completely absorbed into the r\ m factors. On the other hand, the temperature 
dependence of $ and E is set by Ti ns and is weak at T ~ A s . Estimating 77 using E 2 (w) + 
$ 2 (uj) ~ A 2 , we find 77 — 1 ~ T/A s . We see that at low T <C Ti ns the system behaves as a dirty 
superconductor, the role of 7/ A ratio ( where 7 is the elastic scattering rate due to impurities) 
is played by T/A s . Using this analogy and the results for dirty superconductors [[l6],[l7j we find 
that the density of states and the two-particle response functions at finite momentum, such 
as II (u>), are unaffected by T/A s ratio as they have the same form in terms of $ m and E m 
as in terms of 4> m and E TO . On the other hand, the single particle spectral function and the 
two particle response functions at zero external momentum (i.e., the Mcissncr kernel and the 
superfluid stiffness D s ) scale as -q -1 and are substantially reduced above T ~ A s . As in dirty 
superconductors, this reduction can be absorbed into the renormalization of the quasiparticle 
mass m* /m = dE/ 'du> ~ 77 and is not associated with the reduction of the superfluid density. 

These features are present in our numerical solution of the Eliashberg equations. In Fig. [I] 
we show representative results for Re$(ti;) and ImE(aj) for different T. We clearly see that 
at T ~ A s , sharp structures at lo = A + A s transform into broader structures at w = A. 
The latter are due to the fact that in real frequencies, rj(u>) is peaked at u> = A, and the 
amplitude of this peak increases with T In Fig. || we show the behavior of the dynamical spin 
susceptibility and the fcrmionic spectral function (multiplied by np(io)) for different T. The 
insets show the superfluid stiffness, D S (T), extracted from the computation of the Meissner 
kernel, and the NMR relaxation rate \jT\T cx n(w)/a;| w _ >0 . We see that the residue of the 
peak in the spectral function, and the superfluid stiffness sharply decrease above A s . On the 
other hand, 1/TjT and the peak in x'q ( w ) are much less sensitive to the ratio T/A s . We 
explicitly verified by analyzing larger A = 20 that \/T\T does not change much at T ~ A s . 

To this end, we therefore find a crossover in the system behavior at T ~ A s , similar to a 
crossover at 7 ~ A in dirty superconductors. The superfluid stiffness is reduced at T > A s 
due to mass renormalization. However, just as in dirty superconductors, this reduction of the 
stiffness does not give rise to a substantial reduction to T c due to phase fluctuations simply 
because still D s ~ EF(m/m*) 3> A,T C . This follows from the fact that within Eliashberg 
theory E(w - A) - r]A < E F f§. 

We now argue that another physics emerges at T < Ti ns and yields an extra reduction 
of D s unrelated to a mass renormalization. Indeed, we clearly see in Fig [j] that the peak 
in Re<I>(a;) disappears at T — To < T jns . Between To and Ti ns , Re$(o;) ^ 0, but it mono- 
tonically decreases with frequency. Simultaneously, ImE(w) becomes roughly the same as in 
the normal state, the peaks in the dynamical spin susceptibility and in the fermionic spectral 
function become very broad, X/T\T nearly reaches its normal state value, and the superfluid 
stiffness virtually disappears (see the insets in Fig. ||). This behavior is not associated with the 
impurity- like scattering of thermal fluctuations as then the peak in Re$(w) at oj — A would 
disappear only at T,- ns . Rather it implies that immediately below Ti ns there is no feedback 
from pairing on fermionic excitations in the sense that strong ZroS(w = 0) in the normal 
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Fig. 3 - a). The phase diagram emerging from the solution of the Eliashberg equations. Without 
phase fluctuations, T = T lnB is a transition line, while the other two are the crossover lines. Above 
A s , superfluid stiffness is reduced due to enhancement of the effective mass, like in dirty supercon- 
ductors. Above To, it is further reduced due to a reduction of the superfluid density. We conject that 
fluctuations destroy superconductivity at T ~ To < T ins . b) The cuprate phase diagram that follows 
from Fig. a). 



state does not convert into a mass renormalization at T < Ti ns ■ The reduction of D s without 
corresponding mass enhancement likely implies that the superfluid density n s {T) is reduced 
compared to what it would be for a BCS d— wave superconductor at this temperature. We 
conject that due to this extra strong reduction of D s , fluctuations of the gap destroy coherent 
fermionic pairing at T ~ T < T ins [Q. Still, however, the leading edge gap in the spectral 
function remains finite for all T < Ti ns and just fills in by fluctuations at T ~ To. 

The phase diagram which emerges from our consideration is presented in Fig. [| Without 
fluctuations, there is a true transition at T = Ti ns , and two crossover lines at T ~ To and 
T ~ A s . For A > 1, A s < T < T ins . We found that T and A s merge at A - 2, but still both 
remain smaller than T; ns . For smaller A, the distance between T and T[ ns gradually decreases 
and eventually disappears. With fluctuations, it is likely that coherent superconductivity 
appears only at T ~ T . We cannot argue definitely whether this implies that T ins becomes a 
crossover temperature below which the system begins creating disordered singlet pairs which 
condense at T , or there is an Ising-like transition at T ~ T ins where singlets are ordered 
into columnar dimers. The first possibility is, in our opinion, a realization of Anderson's 
RVB ide a fl9|| . The second possibility is a realization of the scenario suggested by Vojta and 
Sachdev [|20f . 

We conclude the paper by summarizing what we obtained. Our key result is the phase 
diagram in Fig.||. We found that there are two different fluctuation effects which govern the 
system behavior below the onset of the pairing instability. First effect is due to thermal spin 
fluctuations, and its role is equivalent to that of nonmagnetic impurities in s— wave super- 
conductors. This effect may account for the reduction of the superfluid stiffness, but cannot 
destroy superconductivity. The second effect is a quantum one - we found that below T iT j S , 
there is no immediate feedback effect on fermionic self-energy, and low-frequency fermionic 
excitations remain overdamped despite pairing. In this situation, singlet pairs still diffuse 
rather than propagate, and superconducting condensate is easily destroyed by fluctuations. 

The data for cuprates seem to indicate that T c scales with the superfluid stiffness at 
T = pl[ and also with the resonance neutron frequency A s |l3| . To account for these results 
in our theory, it is necessary that T and A s coincide. This does happen at intermediate A 
(and our theory does predict that near optimal doping T c scales with A s ) but not at A 3> 1. In 
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other words, the present theory underestimates quantum fluctuations at very strong coupling. 

A final remark. Our results bear some similarities but also some discrepancies with the 
results of the Eliashberg study of phonon superconductors at vanishing Debye frequency p2| . 
The comparison with phonon case requires a separate study and will be presented elsewhere. 
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